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Abstract 

Our velocity measurements on the quasi-two-dimensional turbulent flow in 
a rapidly rotating annulus yield an inverse energy cascade with E{k) ~ 
rather than the expected E{k) ~ k~^^^. The probability distribution func- 
tions for longitudinal velocity differences, 5v{t) = v{x + t) — v{x), are 
self-similar (scale independent) but strongly non-Gaussian, which suggests 
that the coherent vortices play a significant role. The structure functions, 
{[Sv{£)]P) ~ i'^p, exhibit anomalous scaling: Cp = f r^-ther than the expected 
— 3- 
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In large scale flows in the earth's atmosphere and oceans or in gaseous planets, the 
Coriolis force dominates other forces and to lowest order is balanced by pressure gradi- 
ents (geostrophic balance). The dimensionless number that characterizes this regime is the 
Rossby number, the ratio of magnitudes of the inertial term in the Navier-Stokes equation 
{u- Vm, where u is the velocity) to the Coriolis term {2Q x u, where Q is the rotation rate). 
For planetary flows on large scales the Rossby number is typically in the range 0.05-2, and 
turbulence in such flows is quite different from three-dimensional (3D) turbulence in an iner- 
tial reference frame. Despite the importance of planetary flows, geostrophic turbulence has 
scarcely been examined in laboratory experiments. Our experiments in a rotating annulus 
are the flrst to determine the scaling properties of turbulence in a low Rossby number flow. 

One of the most signiflcant effects of rapid rotation on a fluid is the two- 
dimensionalization of the flow [|I|. Recent numerical simulations of a rotating turbulent 
flow recovered quasi-2-dimensionality even in the case of 3D forcing [0 . Turbulence in a 2D 
rapidly rotating flow has some similarities to other 2D turbulent flows, but the rapid rotation 
leads to the formation of large vortices and jets, hence a flow that is more inhomogeneous 
and anisotropic 

The statistics of velocity differences for 2D turbulence are predicted to be self-similar , 
i.e., the probability distribution function (PDF) for the difference between velocities mea- 
sured at two points, 6v{i) = v{x + i) — v{x) (where v is along the line connecting 
the two points), should have a functional form independent of the separation i. Self- 
similarity in a turbulent flow can also be determined from the scaling of structure functions, 
Sp{i) =< [6v{i)Y >~ for any self-similar flow, (p will vary linearly with p [||. In partic- 
ular, Kolmogorov's 1941 theory for homogeneous isotropic turbulence predicts the existence 
of an inertial range with = |. 

Experiments and numerical simulations have shown that 3D turbulence is not self- 
similar |jl|. The observed deviations from a linear dependence of (p on p are attributed 
to the stretching and folding of vortices. As a vortex stretches, it shrinks in lateral extent 
until it collapses onto a singular line. Consequently, the energy dissipation rate varies with 
length scale. Vortex stretching is not allowed in 2D turbulence, and recent numerical simula- 
tions of 2D turbulence yielded self-similar behavior (see also [0,3). A recent experiment 
yielded self-similar behavior with a Gaussian PDF ^ while two other experiments showed 
scale-dependent statistics as will be discussed later. 

Our experiments on turbulent flow in a rapidly rotating annulus yield a self-similar PDF 
for the velocity differences, but the PDF is not Gaussian and the structure function exponent 
values (p differ from the expected p/3. We now describe our experiments and then present 
the results. 

Our apparatus consists of an annular tank flUed with water and covered with a solid lid; 



the inner diameter of the tank is 21.6 cm and the outer diameter is 86.4 cm |I0|]. The depth 
of the tank increases from 17.1 cm at the inner radius to 20.3 cm at the outer radius ||11|| . 
Flow in the annulus is produced by continuously pumping water in a closed circuit through 
two concentric rings of 120 holes each in the bottom of the tank; the source ring is at a radius 
of 18.9 cm and the sink ring is at 35.1 cm. The radially outward flux from the pumping 
couples with the Coriolis force to generate a counter-rotating jet, which is wide and turbulent 
over a wide range of parameters. The rapid rotation leads to a flow that is essentially 2D 
except in the thin Ekman boundary layers at the top and bottom surfaces [pUj]. 
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In the present experiments the tank rotates at 11.0 rad/s, sufficiently fast to produce an 
essentially 2D flow, and the flux is 150 cm^/s, sufficiently large to produce an inverse energy 
cascade. The azimuthal velocity is measured using two hot fllm probes which are inserted 
through the top lid and extend 1 cm into the water on opposite sides of the tank, midway 
between the inner and outer walls. Each probe was sampled at 150 Hz for periods of two 
hours, giving 10® data points per probe for each run, and the measurements were repeated 
four times, yielding a total of 8 x 10® data points. Using the maximum velocity {Umax = 
22 cm/s) as the velocity scale and the distance between the forcing rings (L = 16.2 cm) as 
the integral length scale, we calculate the Reynolds number to be 35, 000 and the Rossby 
number (Umaxf^^L) to be 0.06. 

Instantaneous 2D velocity flelds were obtained using a Particle Image Velocimetry (PIV) 
system with a horizontal light sheet at mid-fluid depth and a rotating camera above the 
tank. For each flow condition, 50 instantaneous velocity flelds were obtained, equivalent to 
approximately 2 x 10^ velocity values at the radius of the hot fllm probes. Though this 
sample size is inadequate for high order statistics, the spatial information provided by the 
PIV measurements complements the long velocity time series obtained with the hot fllm 
probes. We also made PIV measurements at the same rotation rate at a higher pumping 
rate of 550 cmVs {Re = 100, 000). 

Two vorticity proflles obtained with the PIV system are shown in Fig. |I[ When the 
pump is flrst turned on, the flow consists of rings of cyclonic and anti-cyclonic vortices that 
form above the outlet and inlet rings, respectively. Vortices of like sign immediately begin 
merging and growing in size, and a counter-rotating jet forms between the two rings. At long 
times the inverse energy cascade leads to large vortices which are larger at higher pumping 
rates, as Fig. |I| illustrates. 

We compute energy power spectra from the velocity time series data assuming Taylor's 
frozen turbulence hypothesis, which is applicable because the turbulent intensity (ratio of 
the rms velocity fluctuation to the mean velocity) is less than 10%. The spectra contain a 
region with E{k) ~ k~'^ (Fig- 0); in contrast with Kraichnan's prediction of E{k) ~ k~^^^ for 
the inverse cascade Energy spectra obtained from PIV measurements of the azimuthal 
velocity data at Re = 100, 000 show a scaling consistent with that obtained from the time 
series data, as shown in the upper curve of Fig. 0. A prediction of E{k) ~ has been made 
for a turbulent flow made up of Lundgren spiral vortices that were not allowed to stretch in 



the third dimension |T2[. Such scaling has also been predicted for rotating turbulent flows 



as they become 2D 

At high wavenumbers, our power spectra are consistent with those previously found for 
the forward (enstrophy) cascade, E{k) ~ where 3 < n < 4. However, our spectral 
range is too small to deduce a value for n. 

By plotting the PDF for the velocity differences 6v{i) (Fig. §), we obtain the flrst indi- 
cation of self-similarity: data for different separations i collapse onto a single curve when 
normalized. This curve is far from Gaussian, and the enhanced probability in the tails is 
likely due to the strong velocity differences that arise as coherent vortices pass the probes, 
as noted by She et al. |l^ . 

Another indication of anomalous behavior is the scaling of the standard deviation of the 
velocity differences, Svrms- Our data suggest Svrms{(^) ~ ^^''^ (see insets of Fig. H) rather than 
the expected 5vrms{'^) ~ the observed scaling of 5vrms is consistent with the scaling of 
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the energy spectrum in the inverse cascade range (Fig. 

Now we examine the structure function scahng, plotting Sp as a function of i, as shown 
in Fig. ^(a). There is a scahng region (labeled A) from about 2 cm to 8 cm, and some 
indication of another scaling region (labeled B) for i < 2 cm; the transition between the two 
regions can be seen more clearly in the plot of Sio/i^'^ vs. i in the inset. The 2 cm lower limit 
of scaling of the structure functions is below the 5 cm {k = 1.2 cm~^ in Fig. H) length where 
the energy spectrum behavior changes; this difference might be due to systematic error in 
calibration of the probes (note in Fig. 2 the difference in two probes at high wavenumbers), 
or perhaps because the scaling regions in real space and in Fourier space do not exactly 
correspond, as discussed in Ref. [|I|] (p. 62). 

One test for the existence of an inverse energy cascade in 2D turbulence is the sign 
of Ss{£): Kolmogorov's four- fifths law can be written for anisotropic flows as e = — |V£ ■ 
(|5v(£) p5v(£)), where v is the full 3D velocity vector and e is the mean rate of energy 
transfer (Ref. [|I|, p. 88). For the 3D forward cascade e is positive; a negative value of e 
corresponds to an inverse cascade. If the anisotropy is not too strong, the longitudinal S3 
dominates the transverse S^, and one can obtain the sign of e from measurements of the 
longitudinal structure function alone. While our flows are anisotropic, we find that PDFs 
for velocity differences from PIV measurements in the radial and azimuthal direction are in 
agreement, indicating that the anisotropy does not strongly affect the velocity differences. 
The hot film data yield 5*3 > for i < 10 cm, which suggests that the inverse energy cascade 
stops at that point. This length corresponds visually to the scale of the largest vortices; thus 
region A corresponds to the inverse cascade. Scales below about 1 cm (the distance between 
the holes in the inner ring) presumably correspond to the forward enstrophy cascade. 

The now standard way of extracting structure function exponents from data with limited 
inertial range is the technique introduced by Benzi et al. called Extended Self-Similarity 
(ESS), where the (p values are obtained from the slope of log-log plots of Sp vs. S3. The 
ESS plots of our data. Fig. |^(b), exhibit power law scaling through both ranges A and B, 
for lengths 0.5 < i < 15 cm. Since 5*3 should contain the same trends as the other structure 
functions, it is not surprising that the switching from region A to B is lost by this method. 

The exponent values (p deduced from the plot of Sp{i) [region A in Fig. |(a)] and from 
the ESS plot [Fig. ^(h)] are compared in Fig. |[ Region A of Fig. ^(a) yields Cp,A = p/2, in 
contrast with theory, although Cp,Ess = p/3, as it must for a self-similar flow. 

Two recent experiments on turbulence in a quasi-2D fluid, a soap film, found departures 
from linear scaling of (p vs. p [|^Q , indicating that these flows were not self-similar. Another 
experiment on a quasi-2D flow, a magnetically driven shallow layer of electrolyte, yielded an 
inverse energy cascade with self-similar behavior as indicated by both the PDFs for 6v and 
by the scaling of the structure function computed using ESS, which scaled as (p = p/3 0; 
however, the PDFs for 6v were close to Gaussian, in contrast to those yielded by our flow 
(Fig. ^). A recent numerical simulation of 2D turbulence also yielded self-similar PDFs for 
6v, where the departure from Gaussianity, though small, contained significant information 
about the flow @]. 

The p/2 scaling obtained for (p in our flow may be a consequence of an inverse cascade 
driven by a radial velocity shear. This process would yield an energy flux e ~ SvsiSvw)"^ / i, 
where Vg would be a radial velocity due to the interaction of the vortices with the mean 
shear, v\\ the azimuthal velocity and i the azimuthal separation. For a scale-invariant e, the 
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resulting p*'* order velocity structure function would have the scaling exponent p/2. 

In summary, we find that turbulence in our rapidly rotating flow is self-similar, as in- 
dicated by the collapse of the PDFs onto a single curve and by the linear scaling of the 
structure function exponents. This self similarity extends over both the inverse and forward 
cascades. In the inverse cascade range, the flow exhibits anomalous scaling of the energy, 
E{k) ~ A;^^, and the structure functions, Sp ~ ^^Z^, but this anomalous scaling is missed in 
an Extended Self Similarity analysis. We conclude that the rotation and strong shear make 
this flow different from 2D turbulent flows in non-rotating systems. Our results suggest that 
the transport and mixing in low Rossby number flows such as the atmosphere and oceans 
may be different from these processes in isotropic homogeneous 2D turbulent flows. 

This research was supported by a grant from the Office of Naval Research. Z-S. S. 
acknowledges partial support from the Minister of Education in China and the Natural 
Science Foundation of China. 
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FIGURES 




FIG. 1. Vorticity field at Reynolds number 35,000 (Q = 150 cm^/s) and 100,000 (Q = 550 
cm^/s) for rotation rate 11.0 rad/s. Vortices with light (dark) center are cyclonic (anticyclonic) . 
The vortices are advected by the clockwise jet. 
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FIG. 2. Energy spectra (arb. units) with a dotted line showing the Kraichnan k inverse 
cascade and a solid line showing k~'^ behavior. Fits of all of the hot film spectra in the range 
0.1 < A; < 1.22 cm~^ give a slope of —2.04 ± 0.06. The sharp spectral peaks correspond to 
harmonics of the tank's rotation rate, not to dynamics of the flow. 
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FIG. 3. Normalized probability distribution function for the velocity differences, demonstrating 
self-similar behavior: data for different separations {^ = 0.6, 4.6, 9.2, 17.3 cm) collapse onto a 
single curve. The velocity differences are normalized by their standard deviation {Sv)rms and the 
probability by its maximum value, Pmax- The standard deviation of the velocity differences (see 
insets) scales as where a = 0.50 ± 0.06 for Re = 35, 000 and a = 0.54 ± 0.06 for Re = 100, 000. 
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FIG. 4. Even order structure functions (a) as a function of ^ and (b) as a function of 6*3 (an 
Extended Self Similarity plot), for the hot film data at Re = 35,000. The graph of Sm/l^'^ in the 
inset in (a) emphasizes the sharpness of the bend in Sio at £ ~ 2 cm. In (b) the dark symbols 
indicate the region from which the values of C,p were extracted. 
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FIG. 5. Scaling exponents ("p ^-s a function of p for region A of Fig. 0(a), and from the Ex- 
tended Self Similarity analysis. Fig. §(b). Region A, the inverse energy cascade region, yields 
Qp = (0.50 ± 0.03)p, while the ESS scaling shows the expected Cp = (0.33 ± 0.01)p scaling over 
the forward and inverse cascade ranges. The error bars show the standard deviation of the eight 
separate data sets of 10^ points. 
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